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Abstract. The operator approach to analytical evaluation of multi-loop Feyn- 
man diagrams is proposed. We show that the known analytical methods of evalua- 



Qj . tion of massless Feynman integrals, such as the integration by parts method and the 
^ ! method of "uniqueness" (which is based on the star-triangle relation), can be dras- 
I tically simplified by using this operator approach. To demonstrate the advantages 
of the operator method of analytical evaluation of multi-loop Feynman diagrams, 
we calculate ladder diagrams for the massless 0^ theory (analytical results for these 
p . diagrams are expressed in terms of multiple polylogarithms) . It is shown how oper- 
1^ ! ator formalism can be applied to calculation of certain massive Feynman diagrams 
and investigation of Lipatov integrable chain model. 



(N 

a^ ; 1 Introduction 

^ , Elaboration of methods of analytical evaluation of multiple integrals visualized as 



multi-loop Feynman diagrams is motivated from the point of view of physical appli- 



O I cations and the mathematical point of view. 



From the physical point of view such investigations are important, since in cal- 
culations of the physical characteristics in quantum field theories the number of 
Feynman diagrams, in higher order of the perturbation theory, grows so quickly [1] 
rS I that numerical calculations are not sufficient to obtain the desirable precision for 
I the corresponding sums of perturbative integrals. 

From the point of view of mathematics and mathematical physics such investi- 
gations are interesting since they reveal some structures specific for the theory of 
quantum integrable systems [2], [3], [4] (see also [5] and references therein). Besides, 
as a rule, analytical results are expressed in terms of multiple zeta-values and poly- 
logarithms - special functions which are extremely interesting and important objects 
of investigations in modern mathematics and mathematical physics (see, e.g., [6], 
[7], [8]). 

Finally, as we will see in this paper, the analytical results in calculations of 
Feynman diagrams (FDs) give explicit expressions for Green functions of certain 
special quantum mechanical models. And vice versa, the expansion of such Green 
functions over coupling constants gives explicit expressions for multi-loop FDs. This 
connection of perturbative Feynman integrals and Green functions for integrable 
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quantum mechanical models is one of the substantial implications of the proposed 
algebraic approach. 

The main idea of our algebraic method for evaluation of FDs is that we replace 
manipulations with multiple integrals by manipulations with the corresponding alge- 
braic expressions. In other words, identical transformations of multiple perturbative 
integrals are substituted with transformations of generators of special infinite dimen- 
sional algebra. This drastically simphfies all calculations. We stress that here we 
mostly consider the massless case of FDs (see, however. Subsection 4.3). 

2 Feynman diagrams in configuration space 

The Feynman diagrams, which will be considered in this paper, are graphs with 
vertices connected by edges (propagators). To each edge wc assign a complex number 
(the index of the propagator). With each vertex we associate the point in the D- 
dimensional space while the edges of the graph (with index a) are associated 
with the propagator of massless particle 



where x,y & H^, {x — yY" := {J2i=i{xi ~ Hi) {^i — Ui))'^, and a e C. Moreover, 
we consider the graphs with two types of vertices: boldface vertices • denote that 
the corresponding points are integrated over R^. These FDs are called FDs in 
the configuration space. Note that more commonly known the momentum space 
FDs are dual graphs with respect to the graphs which correspond to the FD in the 
configuration space. That is, loops in the momentum space FD are replaced by 
boldface vertices in FD in the configuration space, and vice versa, vertices in the 
momentum space FD are replaced by loops in FD in the configuration space. 

Consider examples of FDs in the configuration space and present the correspond- 
ing multiple integrals. 

1. Graph with 5 vertices and 5 edges (3-point function): 
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2. "Star" graph: 




3. Propagator type graph (2-point function): 
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d'^z d^u d^y d'^w 

(x-z)^«l 2^«2 (z-«)^«3 «^«4 («-2/)^«5 y^"6...{w-xy^°'9 

The problem (which we need to solve when analytically calculate multiple in- 
tegrals corresponding to FDs) consists in searching for special transformations of 
graphs (FDs) such that the number of boldface vertices (the number of integrations) 
decreases at each step. In the next section, we discuss these special transformations 
and describe the corresponding operator formalism, which gives us a possibility to 
represent these transformations using a more compact algebraic language. 




3 Operator formalism 

3.1 Algebraic manipulations with perturbative integrals 

Consider the D-dimensional Euclidean space with the coordinates Xi, {i = 
1,2, . . . , D). We use the notation: x^" — {J^iLi^f)" ■ Let the operators of co- 



ordinates Qi — qj and momenta Pi — pj be generators of the Heisenberg algebra 

[Qk, Pj] = iSkj . 

Introduce the vectors \x) = \{xi}), \k) = \{ki}) that are eigenstates of the oper- 
ators of coordinates and momenta, respectively: qi\x) = Xi \x), Pi\k) = ki\k). We 
normalize these states as follows 

(^1^) = (27r)g/2 ^M^kjXj) , j d^k\k) {k\ = i ^ j d^x\x) {x\ . 
Consider the heat- kernels ("matrix representations") of the operators p~^^: 



where 



13' — D/2 — j3 and T[j3) is the Euler gamma-function. Formula (3.1) relates 



the propagators for massless particles and pseudo-differential operators p "^^ . For 
the operators g^" the "matrix representations" have the form: 



{x\e'^\y)^x''^5''{x-y). (3.2) 



Below we consider three (a,b,c) algebraic relations which are operator analogs 
of relations used for the analytical evaluation of multi-loop perturbative integrals 
for FDs. Recall that these relations give us a possibility to reconstruct FD in such a 
way that the number of integrations (the number of boldface vertices in the graph) 
decreases to zero (this will indicate that a given FD is calculated analytically) . 
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a. Group relation. Consider a convolution product of two propagators: 



/ 



d-z _ G(a',p') /c(„_^)^4^y (3.3) 



We graphically represent this convolution product as 

X , = G{a',P')- x^liz3_y 

Thus, relation (3.3) describes such a reconstruction of the graph in which the number 
of integrations (boldface vertices) decreases by one. This relation is the "matrix 
representation" of the operator identity (group relation) 

p-2a' p~2l3' ^ ^-2(a'+^'). (3.4) 

Indeed, by using (3.1) and (3.2) we easily demonstrate that the "matrix" analog of 
(3.4) 

d^'z \z) {z\p-^^'\y) = 



/ 



coincides with relation (3.3). Note that in the operator relation (3.4) the tedious 
coefficient G{a', (3') is vanished. 

b. Star-triangle relation. This relation is in the basis of the so-called "method 
of uniqueness" [10] (see also [11]) - an efficient method of analytical evaluation of 
FD. In fact, this relation is a special case of the Yang-Baxter equation [2], [13]. The 
star-triangle relation (STR) has the form 



{x - zf-' {z - yyP' {^xfP - y)2(f (y)2a ' 

and was initially used in the framework of investigations of multi-dimensional con- 
formal field theories [12]. The identity (3.5) can be graphically represented as 






= G(a,/?)- 

X (a + /3)' y 

Thus, STR (3.5) describes such a reconstruction of the graph for which the number 
of integrations (boldface vertices) decreases by one. The operator version of this 
relation was proposed in [14] and is written in the form 

Again we note the absence of the coefficient G{a, (3) . To demonstrate the equivalence 
of (3.5) and (3.6) we act on (3.6) by the vectors {x\ and \y) from the left and right, 
respectively, and use the representations (3.1), (3.2). Identity (3.6) can be written 
in the form of the Yang-Baxter equation: 

Ri2{oi) R23{a + (5) RuW) = R23W) Ri2{a + (5) R2z{a) , 
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where Rab{a) = - p^^))^" and [gf \ ] = i % (i, j = 1, . . . , D). Now a few 
remarks about STR (3.6) are in order. 

Remcirk 1. The algebraic version of STR is equivalent to the commutativity for 
the infinite set of operators H{a) —■p^°'q^°': 

H{a) H-\-/3) = H-\-/3) H{a) ^ p^^q^i'^+P) p^^ = p2(a+/3) ^2a _ 

Remark 2. Here we present the algebraic proof of STR (3.6). Introduce an inversion 
operator TZ^ which obeys the conditions 



Xi 



>(A-D) 



TZ^^l, Tl/^qiTlj^^ qi / q , {xi\'R/^^ {—\x 2 



x^ 



ni = n^e^ , Kf^ := TeAP,- 7^A = fPj - 2 (gp) + i{D - A)qj , 

7^A p2/^ 7^A = g'(^+^) p'^ . 



(3.7) 



(3.8) 



We note that TZa — TZq ^ ^ , where TZ = T^-q. Using (3.7), (3.8) the algebraic version 
of STR is proved as follows: 

Remcirk 3. Prom operator identity (3.6) one can immediately obtain new STR 



(p2° + ag-2")(g2(a+/3) + 5p-2(a+/3) ) (^2/3 ^ ^^-2/3) _ 
= if" + cp-2/3)(p2(a+/3) + 5 ^-2(a+/3)^ (^2a + Q^-2a^ ^ 



(3.9) 



where a, 6, c are arbitrary constants. 

Remark 4. One can introduce one more "local" STR [15] which is related to the 
a-representation of perturbative Feynman integrals 

W{f\ai) W{p'\ —) W{f\as) = W{f\ ^) W{f\P2) W{f\ h , 
a2 P3 Pi 

where W^(x^|q;) = exp (— a;^/(2Q!)), and the parameters and are related by 



the identity a, 

resistances in electric networks. 



^ife+/3i/33+fefe 

f3^ 



which is known as star-triangle transformation for 



c. Integration by parts rule [9]. 

First, we present the graphical version of this rule 






Fig. 1 
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With the help of this rule we obtain the reconstruction of graphs in which the number 
of integrations (boldface vertices) does not decrease. However, this rule is extremely 
useful, since the corresponding reconstruction of the graphs leads to variations of 
the indices on the lines, which further permits one to apply previous relations a,b 

and decrease the number of integrations. 

The operator version of the integration by parts rule (Fig. 1) has the form 

[(52 ^2(a+l)l r-2 ^2(/3+l)l 

(27-a-/3)p-gV^ = 4^^+ 1/ 4^^+ 1/ (3-10) 

where a = —a[, 7 = —0:2 and (3 = —a'^. Identity (3.10) can be directly proved by 
using the relations for the Heisenberg algebra: 

[^2^^2(a+l)]^4(^+l)(^ + ^)^2a^ 

where H := ^{piQi + (jiPi) is the dilatation operator. It follows from (3.11) that the 
operators {q^,p^,H} generate the algebra sl{2): 

[f,p']^4H, [H,f]^2f, Hp'^-2p'. (3.12) 
3.2 Group 7i and characters on 7i. 

Note that all operator identities a,b,c (3.4), (3.6), (3.10) are relations on the op- 
erators p^", q'^^ and their products. So it is natural to introduce a group 7i which 
is generated by the operators (Va,/3 G C). The dilatation operator H 

belongs to the group algebra of H in view of the first relation in (3.12). Consider 
any element of the group Ti. 

^(a.) = p-2°'i g-2"2 p-2a', ^-2a, p-2a', . . . ^-2a,, p-2a'^,+, _ (g^^g) 

For a;2m = I^m=o ^2m+i ^ SMch elements form a commutative 

subgroup Ho in 7i (see Remark 1 in Sect. 3). The element (3.13) can be interpreted 
as an operator version of the 3-point function: 





(2;|*(aOlz/) 




2fc-i 2^ "2/s+i y 



Indeed, the corresponding multiple integral is obtained from the representation 

{x\'^{ai)\y) = {x\fr^"''^ q-^""' p-"^"'^ q-^"^ p-2"5 . . . g-2«2fe p-2<+i 

it it it (3.14) 

Jd\ I z, ) {z, I Jd\ I z^ ) {z^ I /ci\ I z^ ) {z^ \ 

if we take into account eqs. (3.1), (3.2). Note that the expression 
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where Xoicti) is a coefficient function, gives the representation for the 2-point func- 
tion, or for the propagator-type FD. The advantage of the operator approach to 
the evaluation of the 3-point function (3.14) consists in that we can apply relations 
a,b,c (3.4), (3.6), (3.10) to the element (3.13) directly at the operator level, i.e. 
use the commutation relations in the group algebra of H instead of making the 
corresponding manipulations with multiple integrals. 

A remarkable fact is that one can define a trace for the elements of the group H 
(generated by the operators {p^", q^^})- First, recall that the dimension regulariza- 
tion scheme requires the identity [17] 



/ 



X 



= Va^O. (3.16) 



The extension of the definition for the integral (3.16) at the point a = was proposed 
in [18] and has the form 

d X 



/ 



X 



D/2 iarg(a) 

where = p^^^^ is the area of the unit hypersphere in R-^, a = \a\e'"''^°' and 6{.) 
is the one-dimensional delta-function. Consider the formal integral for the 2-point 
function (3.15): 

d'^ X d'^x 
-^{x\ ^(a.)|x) = xnia,) J = 

= Jd'^xixl ^{ai)q-^'^\x) =: Tr (^(aOr^^) nQJUP + 7 - D/2\) . (3.18) 

Here P = J2i cti — k^ and Xoicni) is a coefficient function for propagator type diagram 
which now can be interpreted as a character for the group element 

p q ■■■p q ^ eUo , 

where /3 7 = f (or Em=i ^^^m + 7 = Em=o '^2m+i)- The cychc property of the 
trace (3.18) Tr{AB) — Tr{BA) can be checked directly. The definition of the trace 
(3.18) permits one to reduce the evaluation of propagator- type FDs (and searching 
for their symmetries) to the evaluation (and searching symmetries) of vacuum FDs 
(for details see [18]). 

4 Applications 

4.1 Ladder FDs for (f/^ theory in D = 4; relation to conformal 
quantum mechanics 

Consider dimensionally and analytically regularized massless perturbative integrals 



DL{po,Pj^.„p;a,(3,^) 



L 

n/ 



d Pk 

ik=i-' Pk"iPk-pY^ 



n ' (4-1) 

m=0 \rm+i I'm,) 
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which correspond to FDs {xi = po, X2 = Pl+i, X3 = p) 




XI-X3 


VI- V 


Pl -p 


X2-X3 


1 — 


P21 






Pl+il 


X2 








Pi 




X2 



(Pmk = Prr. - Pk) Fig.2 

Perturbative integral (4.1) can be graphically represented in two ways - as diagrams 
in configuration and momentum spaces, as it is shown in Fig. 2, where 7 are 
indices on the lines in the left diagram (in configuration space) , while in the right 
diagram (in momentum space) the indices Pi indicate momenta flowing over the 
lines. As we have mentioned above, these diagrams are dual to each other (boldface 
vertices in the left diagram correspond to the loops in the right diagram). The 
operator version of the integral (4.1) follows from the representation for the left 
diagram 

DL{xa;a,P,j) = {a{J'))-'^-'{x^\p-'^' (uLi ?"'"(?- ^3)"'^^') ^2) • 
It is convenient to consider the generating function for the integrals Dl (4.1) 

9 /ail') 



Dgixa, a, P, 7) = 0(7 ') ^ DL{xa] a, P, 7) = {xi \ if^'- 



X2 



L=0 



If the indices on the lines are related by the condition a + (3 — 2^' 



(4.2) 

, then by using 



properties (3.7), (3.8) of the inversion operator TZ we obtain (for details see [14]): 
i5,(x„;27'-/3,/3,7) = 73^(^1 fp'"' - ^1 ' I ^) > (4-3) 



{xixl) 



where g,,^ = j^^^, = - = _ ^ = _ _ _ ^ ^1,^ ease 

when we fix the indices on the lines as 7' = /3, the generating function Dg (4.3) is 
related to the conformal function 

Dgixa, P,P,D/2-P) 



{u\ \p 



2/3 9f3,l3 



2/3 



(4.4) 



and then taking P — 1 we obtain the Green function for £)-dimensional conformal 
mechanics 



DgiXa, 1,1,1^/2-1) 



(4.5) 



= a(l) 5]/Dl(x„;1,1,D/2-1) . 

L=0 

Thus, we have shown that with a special choice of indices on the lines a — (3 — 1, 
7=y — 1 = 1 — e the ladder diagrams (in momentum space): 
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Xi-X3 1 



X2 - X3 




Fig. 3 



are related to the Green function for £)-dimensional conformal mechanics. Moreover, 
according to the definition of the generating function Dg, the expression Dl for the 
ladder diagram with L loops (L boldface vertices for FD in the configuration space) 
is obtained in the expansion of the Green function (4.5) over the coupling constant 
g (coefficient in order g^). 

The operator method of evaluation of Green function (4.5) is based on the re- 
markable identity [14] 



1 



p2 - g/f 



L=0 ^ ^ 



z2a 



{H-1) 



-2a 



{H-1 + q;)^+i p2 



(4.6) 



where we have used the notation [. . .Jq-i- = jj (d^ [■■■]) The leading terms in 
the expansion of (4.6) over g give the identities 



— 2 — 2 --2 

p q p — 



HH-r 



[log{f), p 



-21 



\-3 



p-' q-'p-' r'p-' = ^^^^ [log(^^), [logm, p-']] + J^' [log(g2), p-'] , 

which can be proved directly and immediately lead to analytical expressions for the 
ladder (box) diagrams with one and two loops {L = 1,2). In general, taking into 
account the integral representation for the rational function of H (in the right-hand 
side of (4.6)) 



(H-l] 



{H-1 + a)^+i 



and using obvious properties of the operator e*^: e*^^"*"?^ \x) — \e we can 
rewrite Green function appearing in (4.5) in the form 



{u\ 



E 

L=0 



L! 



(4.7) 



where 



POO 

/ dtt^ 






JO 







[u — e~^vy 



u 



2(13/2-1) 



[uv } 



(4.8) 
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Thus, the result for the evaluation of the L-loop ladder diagram (Fig. 3) is 

D / 1 \ 2(D/2-L-l) 2 / \ 



2 , , „ ,^2 



where = v'^ = {u - = ^ and Xab = Xa - x^. 

For D = 4 - 2e the function *l(4, 2^) (4.8) is expanded over e 



^ , v"^ {UV)\ ^ ^(k), 



where Zi + Z2 — 2{uv)/u^ and Z1Z2 — v'^/u^. The coefficient functions are 

■ti;' = E ''t " ^' i: (-)•" q" (^■. 2l + ; - /) , (4.10) 

/=0 J-K^-J)- m=0 

where C/" = " ^| is a binomial coefficient and 

t <!(/— m)! 

1^1 -^2 j no,...,nm=l (Ei=o^*) 



is expressed in terms of multiple polylogarithms 



^ ^ ,,,,, 

l^imo,rm,...,mAwo,Wu...,Wr)^ }^ -j—^^ — . (4.11) 

no>ni>->nr>0 '^0 'n ■ • ■ 



The first coefficient (for D = 4 or e = 0) has the form [19], [20] 



4°H^i,^2) = E^f^T^^^ \nf{z,Z2)[u,,_^{z,) - U,,_^{Z2) 

and is expressed via the standard polylogarithms Li^niw) — J2'^=i The following 
coefficient was calculated in [14]: ^^l\zi,Z2) — 

1 g n\ \n'^~"{ziZ2) \{nU^^, (zi) - Li„^ {zi, 1) - Li„_, {zi, )) - {zi ^ Z2) 



{z,-Z2)^^ (-l)-(2L-n)!(n-L)! 
where the functions IA^^{wq,Wi) have been defined in (4.11). 

Remark. Using conformal and scaling properties of the Green function (4.3) one 
can deduce the representation (cf. (4.8)): 
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where the conformal symmetry requires ¥'^'^\ui,U2) = *^T'2t-'^)(mi,M2). From the 
representation (4.12) we obtain the identity 



u-^-'- " {u\ \p^^-g ' ^■iii 



,2(D/2-7) 



\v) 



f3-j ■ 



(4.13) 



where u 



X\ X3 ' X2 Xl Xl2 ' X13 X12 



the concise notation 
cross-ratio identities 



— and we have introduced 

u;i u;i2 ' x'13 Xl2 

To prove eq. (4.13) it is only needed to note the 



{v'f xlr^xl {u-vf {u'-v'f xl^xl 



2 1 1 
XiX]2, 



2 2 



Now for both sides of eq. (4.13) we make the inverse transformation with respect 
to that used in passing from eq. (4.2) to eq. (4.3). As a result, we rewrite (4.13) in 
the form 



^2(7-i?/2) 
^3 



(xi I - 



27 

g u 2 



g2(27-/3)(^-_ ^3)2/3 



\X2) 



2{'r-D/2) 

•^12 



{xi\ [p 



27 _ 



27 



^2{2^-f3)(^q_Xu] 



2/3 



where u 



and V — 



~2 — 7 



-. Note that this identity is nothing but the relation 



on the generating functions for the ladder diagrams (4.2). Expanding both the sides 
over the coupling constant g, we obtain D-dimensional identities for the L-loop 
ladder momentum diagrams in the order g^: 



2:1-3:3 /3 13 . 



2(7+7L-f ) ^, L(/3-j) 



U 2 X 



"1 /3 /3... 



^2 (3/3. 



X 



2(7+7L-f ) ^ Lifl-y) 



12 



V 2 X 



where (3, (3 = 2j — P and 7' = D/2 — 7 arc special indices on the lines and xi,X2, X3 
parametrize external momenta. These identities, in the special case D = 4 (e = 0) 
and /3 = 7' = /9 = 1, were obtained in [21] and used there for deriving many 
remarkable relations for various planar FDs. 
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4.2 Application to Lipatov chain model 

It was shown in [3] that a wave function for bound states of gluons at high energies 
has the property of the holomorphic factorization. The Lipatov Hamihonian for 
each of these two holomorphic subsystems has the form: H = Yh^i Hu+i, where 

Hik = Pi ln(pife)pri + \n{pik)Pk'^ + ln{piPk) + 27 = (4.14) 

= 2 In(pifc) + In(piPfc) p~^^ + 27 . (4.15) 

Here 7 = — r'(l) is the Euler constant, pik = Qi — Qk, where qi are complex co- 
ordinates, and momentum operators pi — —%-§^. are complex derivatives. In this 
subsection we show how one can easily demonstrate the equality of two expressions 
(4.14) and (4.15) for the Lipatov pair Hamiltonian by using the operator technique 

discussed above. 

Note that the expression (4.15) (up to the constant 27) appears in the expansion 
over e of the following operator: 

^ife(e) p]^'^iPkrp1k^' = 1 + e (2 ln(p,fc) + p,k HPiPk) P^k) + ■ ■ ■ ■ (4.16) 

Now we use the one-dimensional analog of the operator "star-triangle" identity (3.6) 

P^Pt^^pi = P^P^^^P? ^ PtiPt^%i = pfpt^^P? ■ (4.17) 

Then, we have 

Pik'^ iPiPkY Pik^'^ = Pik'^ PI Pik Pik Pk Pik^'^ ~ Pi ^ Pik Pi~^^ Pk Pik Pk ~ 

+ e (pY^ ln{pik) Pi + Pk^ In(pifc) pk + \n{piPkfj + . . . , 

and this proves the equivalence of expressions (4.14) and (4.15). 

We stress that from the point of view of the integrability of the Lipatov model 
the form of the pair Hamiltonian (4.14) and its relation to the operator (4.16) is not 
accidental. Indeed, the i?-operator (4.16) satisfies the Yang-Baxter equation 

Rii+li^) -Ri+li+2(e + e') Rii+li.^') = -Rj+l 1+2(^0 Rii+l{^ + Ri+li+2i^) ; 

which can be easily proved by using the operator "star-triangle" relation (4.17). 
Then the complete holomorphic Hamiltonian H = -f^ii+i appears in the expan- 
sion over e of the monodromy matrix (in the order e^) 

T(l,2...,n+l){e) = Rl2{e) -R2 3(e) Rsii^) ■ • •-Rnn+i(e) . 

Finally, we should hke to note that recently the solutions of the Yang-Baxter 
equation for principal series of infinite dimensional representations of SL{N) have 
been considered in [22]. The results of [22] generalize the discussion presented in 
this subsection. 
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4.3 Diagrams with massive propagators 

In this subsection we consider an example of the operator approach to analytical 
evaluation of the 1-loop 3-point function with one massive propagator. 

First, we use the automorphism (p^ ^ (f, H ^ —H) of the s/(2)-algebra (3.12) 
to write the first relation in (3.11) as 

[g2/3^ p2| _ 4^ _ ^ + l)g2(/3-l) ^ (4.18) 

and, then, generalize it by introducing the massive parameter m as follows: 

[(^2 + ^2)/3^ ^2] =4p(^H + m^drn^ - p + l){q^ + 771^)^^-'^ . (4.19) 

This identity can be converted into the integral form 

p2 p2 4^ ' 

from which the representation for the 3-point function is deduced 

/ I 1 /'^2 , 2xf/3-i) 1 I \ «(1) f°° 1^ trD/2-n (e~*a^i + "^^)^ - (e*x? + m^)^ 

{xih^iq^ + m^yf' '^^ X2) = / dt e*^^/^ ^ ^2(d/2-i) " 

Ap Jo [e^xi — X2y^ ' ' 

Here the left-hand side is represented in the form of the perturbative integral and 
we obtain the equality 



I ik- 



d^'k 



{k - a;i)2(o/2-i)(A;2 + m^y-^^{k - X2)2W2-i) 

= J_ r dt .*(-D/2-i) (^"*^2 + my - {xy + my 
Aa{l)Jo /3(xie*-X2)2(^/2-i) 

Finally, we consider the limit D —> 4, j3 ^ for this relation and deduce the identity 
d^k _ 2 J /e"*x|-Fm^\ 



J (k — x^)'^(k^ + 'm^)(k — Xo)'^ Jo (xt — e~^Xo)^ ^ 



(A; — Xi)2(A;2 -|- m2)(A; — 2:2)^ Jo (xi — e ^x-i)^ \^x\ + w? j ' 

which is important for physical applications and corresponds to the evaluation of 
the 3-point one-loop FD (in the momentum space) with one massive line: 



Xi — X' 




k — X2 
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Conclusion 



A large enough class of perturbative integrals (which are graphically represented as 
FDs) generates a commutative algebra of functions Ai. For some subset Aio G Ai oi 
generators of this algebra, the explicit analytical expressions in terms of special func- 
tions (of the type of multiple polylogarithms) are known. The problem of analytical 
evaluation of a specific perturbative integral / e At is now reduced to searching for 
the representation of / (e.g., by means of transformations (a,b,c) from Section 3) in 
terms of the elements of Aio. Unfortunately, one cannot always succeed in finding 
such a representation, since the subset Aio which is known nowadays, as well as the 
collection of transformations (of the type a,b,c), do not give the possibility to speak 
about A4o as a complete system in A4. In this case, by using the transformations 
(a,b,c), one can obtain certain relations between the elements 1^ E Ai, which arc 
considered as a set of functional equations for la- Thus, the problem of analytical 
evaluation is reduced to the extraction of an independent set of equations and then 
to searching of their solution. 

Now let us make a few remarks about the results presented above. 

1. It should be noted that the coefficient functions ^l{u,v) (4.7) appear in the 
calculations of the 4-point functions in the N — A supersymmetric Yang-Mills 
theory [16]. 

2. The proposed operator relations (3.6), (4.17) not only clarify the structure of 
the integrable Lipatov model (see subsection 4.2), but also help to investigate 
certain generalizations of this model (see [22]). 

3. The important problem is the search for generalizations of the described al- 
gebraic formalism in the cases of supersymmetric quantum mechanics and for 
massive propagators. In the last case we have succeeded in calculating the spe- 
cial 3-point FD with one massive propagator (see subsection 4.3). However, 
this calculation is particular. From this point of view it would be important to 
calculate the coefficients $l(m, "w;?^^) in the expansion over g of the spectral 
Green function for conformal mechanics: 

-j^ oo 

{p^ - g/r + rn^) t^o 

I am grateful to G. Arutyunov, S.E. Derkachov, L.N. Lipatov, O.V. Ogievet- 

sky and E. Sokatchcv for helpful discussions and comments. I would like also to 
thank the Organizers of the Conference " Quarks-2006" (Repino, May 2006) where 
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